Subtraction of the spurious translational mode from the RPA response function 
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It is well known that within self-consistent Random Phase Approximation (RPA) on top of 
Hartree-Fock (HF), the translational symmetry should be restored. Due to approximations at 
the level of the practical implementation, this restoration may be only partial. As a result, one 
has spurious contributions in the physical quantities that are extracted from RPA. While there are 
several recipes in the literature to overcome this drawback in order to produce transition densites 
or strength functions that are free from spurious contamination, there is no formalism associated 
with the full RPA response function. We present such formalism in this paper. Our goal is to 
avoid spurious contamination when the response function is used in many-body frameworks like the 
particle-vibration coupling theory. 

PACS numbers: 21.60.Jz, 21.10.Re, 24.30.Cz 



I. INTRODUCTION 

RPA is one of the basic tools to study the nuclear vi- 
brational excitations. The basic theory is well known, 
and several versions have been implemented since the 
early days of nuclear structure studies. However, while 
in the past the underlying mean field and the residual 
particle-hole (p-h) interactions were often purely empir- 
ical, nowadays self-consistent HF plus RPA calculations 
are feasible with quite sophisticated effective forces. The 
importance of these calculations stems from the possi- 
bility of testing their predictive power, e.g., when going 
towards neutron-rich (or neutron-deficient) isotopes that 
can be studied at present or future radioactive beam fa- 
cilities. 

One of the advantages of RPA is that it is a conserv- 
ing theory, that is, the symmetries of the Hamiltonian 
that are broken at the HF level are in principle restored 
in RPA. The detailed way in which this restoration is 
realized has been already widely discussed in the liter- 
ature. The reader can consult standard textbooks (cf. 
Refs. or review papers like Ref. 0. The purpose 

of the present work is not to re-discuss basic questions. 
We shall focus instead on practical issues, in the case of 
the restoration of translational symmetry (which is asso- 
ciated with the appearance of a zero-energy translational 
mode in the RPA J n = 1~ spectrum). 

In practice, it is impossible to obtain this translational 
mode at exactly zero energy due to some theoretical or 
numerical approximation. If we restrict ourselves to the 
case of nonrelativistic calculations performed with effec- 
tive Skyrme forces, no fully self-consistent calculations 
that treat the continuum exactly are available. Fully 
self-consistent calculations in which the continuum is dis- 
cretized do exist 0, H| , and in this case the energy of the 



spurious translational mode goes to zero as the dimension 
of the box in which the system is set increases. Both in 
this case, and in the case of continuum calculations which 
are the object of our study, one has then to deal with the 
spurious contamination of several quantities. 

In fact, if the spurious state from a numerical calcu- 
lation does not lie at zero energy, and its wavefunction 
does not coincide with the true translational wavefunc- 
tion, the other (physical) states contain some admixture 
of the ideal spurious state. These admixtures should be 
eliminated. Various procedures have been proposed, de- 
pending on the specific RPA implementation. In the case 
of coordinate-space RPA the basic question is how to 
remove spurious admixtures from the strength function 
associated with a given operator F(r) acting in the 1~ 
channel, namely from 



S(E) = J2\(n\F\6)\ 2 S(E-E n ), 



(1) 



where we suppose that |n) are the RPA states having re- 
spectively energy E n and |0) is the nuclear ground-state. 
We assume that the definition can be formally extended 
to the continuum case, when the spectrum is not dis- 
crete. Sometimes one is willing to remove the spurious 
admixtures directly from the transition densities, that 
are defined as 



6p n (r) = (n\p(v)\0), 



(2) 
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where p is the density operator. Once more, this defini- 
tion holds if the spectrum is continuum and the transition 
density is associated with a state at every positive energy 
lo: we will use the notation Sp(ruj) in what follows. 

These two problems have been dealt with extensively 
in the literature (see Refs. [1,0-03 )■ The solutions in- 
clude projection methods, use of corrected effective op- 
erators etc. However, the present work is motivated by 
the use of the dipole RPA spectrum not by itself, but as 
an in put to particle-vibration coupling (PVC) calcula- 
tions [PZHlij . In PVC calculations performed in coordi- 
nate space with proper continuum coupling [l5j the full 
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RPA response function enters; consequently, one must 
deal with the subtraction of spurious contamination from 
that quantity. Since this has not been discussed in the 
literature previously, to our best knowledge, we discuss 
the problem in the present paper. The next Section deals 
with the theoretical formalism. 



II. METHOD 

As recalled in the Introduction, the spurious dipole 
mode is originated by the restoration of the translational 
symmetry in RPA. Translations along one direction, say 
the z-axis, are generated by the total momentum oper- 
ator P which is a conserved quantity, i.e., [H, P] = 0. 
As discussed in 0, 0, [l6[ , the ideal spurious state must 
have X and Y amplitudes associated with the p-h matrix 
elements of P. Because of the symmetry of the RPA ma- 
trix, its non-zero eigenvalues and eigenvectors appear in 
pairs. The partner of the spurious state has amplitudes 
related to the matrix element of the operator Q that sat- 
isfies [H, Q] — -fi^P- This latter equation is associated 
with the Galilean invariance as it is shown in Refs. [HIH- 

In this Section, we show a method which removes the 
spurious contribution from the RPA response function in 
coordinate space representation. We express the formu- 
las that have been introduced in Ref. [17| in coordinate 
space representation, and extend them to case of the RPA 
response function. 

In order to remove spurious admixtures, first of all, as 
in Ref. [l7|, we suppose the validity of 



Sp(uj) = Sp(us) + \p(bj)8p P + Xq(uj)Spq, 



(3) 



together with the orthogonality conditions [6p(uj), Spp] — 
[S p(oj) , S pq] = 0, for the "physical" transition density 
operator Sp{ui) in its second-quantized form. Sp{uj) is 
the second-quantized form of the RPA transition density 
operator. Spp and SpQ are related with P and Q by 

Sp P = *E / dri>l(ra)(-ihV z )i> q (ra)=iP, (4) 
q,a J 

Spq = / dr4>l(ra)zi> q (ra) = iQ. (5) 

q,a J 

The canonicity conditions [8pQ,8pp] = —[Q,P] = —iTiA 
hold. 

The transition density Sp{ru) in coordinate space rep- 
resentation obeys 

^, ja (rw) = (0|[^,, a (r),*p(w)]|0>, (6) 
where p q ^ a {f) is the density operator defined by 
PqA r ) = ^2^l(ra)O a (r)rjj q (ra) 

a 

with 6 a (r) e {l,^± V,...}. (7) 



Using Eq. ([3]) and Eq. ([6]), we obtain 

8p q<a (rw) = Sp q , a (ruj) + \p(uj)(0\{p q , a (r),Spp}\0) 

+\ Q {u>)(0\[pq, a (r),Sp Q }\0). (8) 

By using Eq. (0} and Eq. ((6]) with the orthogonality 
condition and the canonicity condition, \p(w) and Xq(ui) 
are represented by the transition densities in coordinate 
space representation as 



Ap(w) 



Aq(w) 



(0\{Sp Q ,5p(oj)}\0) 
(0\[Sp Q ,Sp P }\0) 

(0\[6pp,6ftu)]\0) 1 



hJ^J drz5 Pi( ru} ) 



(9) 



(0\[$PQ,$Pl 



= lE / dr6 jz,q( rUJ )- 

q J 

(10) 

For example, if O a (r) = 1, then Eq. ([S]) becomes 
Sp q {ruj) = 8p q {ruj) 

+V zPq (r)^J2 J dr'z'5p q ,{r'uj) (11) 
q' 

because {0\[p q (r),Spp]\0) — HV z p q (r). As a further ex- 
ample, if O a (r) , one has 

S2ij (rw) = 52ij (ru) 



+P 9 (r)g [ E / dr'S2ij z , q ,(r'u;) 



(12) 



because {Q\[2ij (r), Spq]\0) — ie z p q (r) (e z is the unit 
vector along z). 

In this paper, we extend the above method to single 
out "physical" quantities labelled by a tilde, to the case 
of the response function instead of the transition density. 
We adopt here the same notation as in Ref. [18]. The 
response function R^, depends on two indices a and (3 
that identify the operators as in Eq. (0. 

As in Ref. [IH , one can write 

*P,H = E/ dr'R^irr'-^fAr'), (13) 
q' 

S2ij q (ru;) = E/ dr'R v qq 7' l {rr';u)f q ,{r'), (14) 
q' 

where f(r) is the external operator. In the previous for- 
mula the superscript 1 refers to the identity operator and 
the superscript V_ to the operator 

We can use Eqs. (fTT|) and (fT3|) : then, for the response 
function we can arrive at 

R q 1 qt (rr';u;) = Rl 1 q ,(rr';oj) 

-V z p 9 (r)^^ J dr"z"R^ q ,{r"r'^). (15) 
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This is our main result. It allows extracting the "physi- 
cal" response function R from the RPA result R. 

In the next Section, we illustrate our method by means 
of some numerical results for the strength function S. In 
fact, by using the "physical" response function defined in 
Eq. (fl~5|) , the "physical" strength function can be written 
as 

qq' 

= S{u) 

-^ Im E JdrP q (r)V zPq {r)/A 



Jdr'J dr"z"R$ q ,(r"r';u;)U(r') 

q'q" 



S{w) - SS{lu), 



(16) 



where S(uj) is the strength function defined by the RPA 
without corrections and SS(uj) is the spurious contribu- 
tion. 

As mentioned in the Introduction, the main point of 
correcting the response function instead of the transi- 
tion densities or transition strength is that the response 
function is an input for more sophisticated many-body 
frameworks like the particle- vibration coupling approach. 
In this case, the static HF potential is generalized to a 
time-dependent (or energy-dependent) self-energy func- 
tion. This self-energy function £ is associated with the 
coupling of single- nucleon states with RPA vibrations; 
consequently, it can be defined in terms of the HF Green's 
function Go and the RPA response function. This defini- 
tion can be found in the references mentioned in the In- 
troduction. However, we have recently tried for the first 
time [la ] to formulate this theory in coordinate space. 
In this case, in order to remove the contribution of the 
spurious state by using Eq. (fT5|) . the self-energy function 
must be expressed as 

E(rer, r'a'; uj) 

— K(r)G (ra, r'a'; u - u}')K(r')iR n (rr'; w'), 
> 2?r 

(17) 

where k is the residual (p-h) interaction. 

Using this self-energy function, we can solve the Dyson 
equation G = Gq + GqY<G in coordinate space, and then 
obtain the perturbed Green's function. The detail of the 
formulas will be in our forthcoming paper |15| . where 
also several results are discussed. In the present work, we 
present some few results about the relevance of projecting 
out the spurious contribution. 



III. NUMERICAL RESULTS 

In this Section, we study the strength functions both 
in the case of the electric dipole and the isoscalar com- 



pression operators. For this study, the nucleus of choice 
is 208 Pb. The equations are solved in a radial mesh that 
extends up to 20 fm. The radial step is 0.2 fm. The states 
up to angular momentum ^=12 are included in the model 
space. The Skyrme set employed is SLy5 [19j. The im- 
plementation of the continuum RPA is exactly as in Ref. 

Also we apply our method to Hartree-Fock plus 
particle-vibration coupling. As an example, we show the 
level density in the case of 24 O. We included the RPA 
electric dipole phonons up to 60 MeV. We discuss the 
spurious state contribution to the level density, extracted 
from the Green's function which is solution of the Dyson 
equation (in the level density we adopt a smoothing pa- 
rameter equal to 0.4 MeV). 



A. El excitations 

In Fig. [TJ we display the strength functions (and the 
associated energy-weighted sum rules, EWSRs) obtained 
by using the following El operators: 



Fei = eJ~]rpYiM(rp), 
P =i 



(18) 



and 



N Z Z N 

F E i = e ^ -^ r nY 1M (r n ) - -r r p Y 1M (r p ). (19) 



The operator defined by Eq. (|T9|) is obtained by remov- 
ing the center of mass from the original operator of Eq. 
(|18p. We apply our method to obtain S(ui), by leaving 
the operator equal to ((T8j) . and we obtain the blue curve 
in Fig. [1] One can compare this result to the El spec- 
trum without the spurious subtraction (red curve), and 
to the one obtained by using the original RPA response 
and the operator (|19l) (green curve). In the lower panel 
of the figure we show the running, or cumulated, values 
of EWSRs. From both panels one can argue that our 
method succeeds to remove the spurious state. In par- 
ticular, the success of our microscopic implementation is 
very clear from the detailed strength of the upper panel, 
as no strength remains in the low-energy region down to 
zero. 



B. IS compression dipole 

Similar accurate results can be obtained if the IS dipole 
compression operator is applied. This operator reads 



Fisi — 



A 

E 



rtY 1M (ri) 



(20) 



without any center-or-mass subtraction. Many authors 
have used, to take care of this subtraction, the hydrody- 
namical ansatz that has been originally proposed in Ref. 
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10 15 20 

Energy [MeV] 

FIG. 1: (Color online) The Skyrme continuum RPA results 
for the El strength of 208 Pb with the SLy5 Skyrme interac- 
tion are displayed in the upper panel. The red curve [labelled 
by (1) in the figure] is calculated by using the El operator 
(|18fl . whereas the green curve (2) is calculated by using the 
El operator with the center-of-mass subtracted [Eq. <|19p ], 
and the blue curve (3) is calculated by using the operator of 
Eq. (|18[) but at the same time with the spurious state elim- 
inated from the response function by using our method. In 
the lower pannel we show the running (or cumulated) values 
of the energy weighted sum rules. The dotted lines are the 
double-commutator values of the energy-weighted sum rules 
for the operators (|18l) and (| ltjf) . The formulas are given in the 
Appendix [Eqs. (|A4[) and QA8[I , respectively] and the values 
are shown in Table [I] 



|. This amounts to using the operator 



Fisi — 



A 

E 



(r>< Y 1M (ri), 



(21) 



the factor 5(r 2 )/3 being often denoted with the symbol 

7]. 

We present our results in a form which is analogous to 
the one of the previous subsection in Fig. [2j The blue 
line corresponds to the choice of the bare operator (I2TJ1) . 
and to the miscroscopic subtraction of the center-of-mass 
motion proposed in our present work. The associated re- 
sult can be compared with the use of the operator (|2"Tj) in 
connection with the original RPA response (green curve) 
and to the RPA result without any subtraction of the 
center-of-mass either in the operator or in the response 
function (red curve). 

The result of the red curve includes the spurious con- 
tributions. Both other curves provide a way to remove 
these contributions. Although similar, they are not iden- 



tical: in this respect, the upper panel of Fig. [2] shows 
to which extent the ansatz (f2~Tj) for the effective operator 
is accurate. The small differences are not visible if one 
looks only at the cumulated EWSRs: both the blue and 
green line converge to the same result. 




20 30 
Energy [MeV] 

FIG. 2: (Color online) Same as Fig[T]for the isoscalar com- 
pression dipole spectrum. The nucleus is 208 Pb and the 
Skyrme set is SLy5. The red curve [labelled by (1) in the 
figure] is calculated by using the operator (|20[) . whereas the 
green curve (2) is calculated by using the operator with the 
center-of-mass subtracted [Eq. (|21[l ]. and the blue curve (3) 
is calculated by using the operator Eq. (|20p but at the same 
time with the spurious state eliminated from the response 
function by using our method. In the lower pannel we show 
the running (or cumulated) values of the energy weighted sum 
rules. The dotted lines are the double-commutator values of 
the energy- weighted sum rules for the operators (|20|l and (I21|) . 
The formulas are given in the Appendix and the values are 
shown in Table [I] 



C. Spurious mode subtraction in PVC 

The level density pij can be defined as 
±1 f°° 

py(w) = — / drImGij(rr,u), (22) 

JO 

where Gij can be either the Hartree-Fock Green's func- 
tion or the perturbed Green's function which is obtained 
by solving the Dyson equation. The spurious state may 
affect the self-energy function that appears in the Dyson 
equation because the definition includes the RPA re- 
sponse function. Therefore, the spurious contribution 
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Operator 


mi 


mi 


Smi 


El (0 


872.3 


1350 


477.7 


cm} 


872.3 


872.3 


0.0 


Compression dip. (I20p 


5.588 x 10 6 


1.368 x 10 7 


8.092 x 10 6 




5.588 x 10 6 


5.588 x 10 6 


0.0 



TABLE I: Values of the energy-weighted sum rule for the El 
and the isoscalar compression dipole excitations in 208 Pb with 
SLy5. These values correspond to the dotted lines in Fig. [1] 
and [2] 5m i is the contribution arising from the center-of- 
mass spurious motion. Definition are given in the Appendix: 
in particular rh\ is defined by Eq. (|A2[) . 



(P 3 / 2 ) 



(Sl/2) °, 



( 1 ) Hartree-Fock 
(2) HF+PVC(0, no spurious 
(3) HF+PVC(0 



24 



O SLY5 




-8 -7 -6 -5 -4 



should be systematically eliminated from the response 
function. 

In the neutron-rich nuclei like 24 0, low-lying El 
strength can be larger than in other systems, and the 
spurious state is not well separated from the other phys- 
ical states, as it is shown in Fig. [3] 




10 20 30 40 50 60 

E [MeV] 

FIG. 3: (Color online) Same as Fig. [Qfor 24 0. The solid red 
curve [labelled by (1) in the figure] is calculated by using the 
operator 1)180. and the dotted blue curve (2) is calculated by 
using the operator Eq. (|18p but at the same time with the 
spurious state eliminated from the response function by using 
our method. 



In FigHl the level density is plotted as a function of 
the single particle energy. The dotted black curve is the 
Hartree-Fock level density. The solid red curve is the 
HF+PVC level density, in which the spurious contribu- 
tion is properly removed. The dashed blue curve is the 
HF+PVC level density, but when the spurious state is 
not removed. We have studied neutron level densities, 
and in particular the upper pannel is the level density in 
the (P3/2) _1 case, and the lower pannel is the (si/ 2 )~ 1 
case. By comparing the solid red curve and the dashed 
blue curve, we can see that the unphysical contribution 
of the spurious state cannot be ignored. 



1 ■ - — 1 -=a 

-25 -20 -15 -10 -5 

s.p. energy [MeV] 

FIG. 4: (Color online) The level density defined by Eq. (J22J) 
in 24 O. The dotted black curve [labelled by (1) in the figure] 
is the Hartree-Fock level density. The solid red curve [labelled 
by (2) in the figure] is the HF+PVC level density when the 
spurious state is removed by applying Eq. (|15[) . The dashed 
blue curve [labelled by (3) in the figure] is also the HF+PVC 
level density, but when the spurious state is not removed. The 
upper panel corresponds to (P3/2) _1 , and the lower panel is 
(si/2) -1 orbit. 



IV. SUMMARY 

While RPA provides, in principle, a theoretical frame- 
work that restores the symmetries which are broken at 
the mean field level, approximations that are often made 
can spoil this feature. This is particularly true for the 
case of the translational symmetry. It has been known 
for several decades that the self-consistency violations, 
the cutoff in the model space and/or other limitations 
of the practical RPA implementations leave the spurious 
center-of-mass motion at finite instead of zero energy and 
give associated contaminations of the wavefunctions and 
physical observables related with the RPA states. 

Many papers have been devoted to this problem, yet 
the subject of this work is the removal of spurious contri- 
butions from the RPA full response function, within the 
continuum RPA on top of Skyrme HF framework. We 
have presented a microscopic formalism that is claimed 
to be able to achieve this removal quite successfully, and 
we have validated this statement by means of accurate 
calculations of the El and isoscalar compression spectra 
of 208 Pb. 

The original motivation of this work is the possibility 
of using the RPA response function for particle-vibration 
coupling calculations. Indeed we have shown, by using 
the example of level densities in 24 0, that the results of 
PVC either with or without the proper subtraction of 
the spurious state can substantially differ. Other many- 
body scheme or reaction calculations require the use of 
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the response function. In all these cases, our proposed 
method can display its usefulness. 



Appendix A: Energy-weighted sum rules of interest 

In general, we define the energy-weighted sum rule as- 
sociated with the strength function S(u)) as 

mi = J dujujS(uj). (Al) 

Consistently with Eq. (fT6|) wc define 

fhi = J dojuiS(uj) = mi — 5mi (A2) 

with 

5m\ — I dujujSS(uj). (A3) 



It is well known that mi can be calculated analitically 
by means of the expectation value of the double commu- 
tator. In this Appendix, we only summarize the useful 
results in the case of Skyrme interactions. 



Therefore, 



rh\ = mi — Smi 



e 2 h 2 9 „ 9 - e 2 h 2 9 Z 2 



-Z +—e 2 bl - 
2m An An 2m An A 

e 2 h 2 9 NZ 9 



2m An A An 



— e 2 bl. 



(A8) 



The mi associated with the operator ([TO)) is given by 

MM 



mi 



e 2 h 2 9 NZ 9 ^ 
2m An A An 



(A9) 



and, as expected, Eq. (|A8| an( l ||A9j) are completely 
equivalent. 

2. The compression dipole operator 



The mi with the operator (|20|) is given by 



= ^_A,4(llr 4 ). 
2m47r x 



(A10) 



Also Smi can be derived easily as 



1. The El operator 

The mi associated with the operator ([T8"l) is given by 
e 2 h 2 9 9 



mi 



Z+ — e z bl, 



2m 47r 47r 



(A4) 



with 



b = \ (ti(l + ^ Xl )+t 2 (l + ^x 2 ) ) . (A:,) 



I = drAur 2 p n {r)p p {r). 



(A6) 



The second term of mi is a correction arising from the 
momentum-dependent terms of the Skyrme interaction. 
5m i is given by 



Smi — 



e 2 h 2 9 Z 2 
2m An A' 



(A7) 



(All) 



and therefore 



rhi = mi — 5mi 

= ^± A {llr i )-^-^-A(5r 2 ) 2 
2m An x 1 2m An x 1 



2m 47r 



A(33(r 4 ) -2h{r 2 f 



(A12) 



On the other hand, the energy- weighted sum rule mi 
associated with (|2"TT) reads 



mi= £i A ( 33(r4) - 25(r2)2 ) 



(A13) 



Once more, Eq. (|A12|) and (|A13|) are completely equiva- 
lent. 
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